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We measure the dimuon charge asymmetry A in pp coUisions at a center of mass energy ^/s = 1960 
GcV. The data was recorded with the DO detector and corresponds to an integrated luminosity of 
approximately 1.0 fb~^. Assuming that the asymmetry A is due to asymmetric B'^ B^ mixing 
and decay, we extract the CP-violation parameter of B° mixing and decay: 

'^^'-B")^^ = ^ = -0.0023 ± 0.0011 (stat) ± 0.0008 (syst). 

l + |eso| 4 
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Ago is the dimuon charge asymmetry from decays of B°B° pairs. The general case, with CP 
violation in both B'^ and B^ systems, is also considered. Finally we obtain the forward-backward 
asymmetry that quantifies the tendency of /i"*" to go in the proton direction and to go in the 
anti-proton direction. The results are consistent with the standard model and constrain new physics. 

PACS numbers: 13.25.Hw; 14.40.Nd 



I. INTRODUCTION 

We measure the dimuon charge asymmetry: 

N++ -N— 



A = 



■N- 



(1) 



in pp collisions at a center of mass energy ^/s = 1960 
GeV. iV++ (N^^) is the number of events with two pos- 
itive (negative) muon candidates passing selection cuts. 
The data was recorded with the DO detector at the Fer- 
milab Tevatron between 2002 and 2005. The exposed 
integrated luminosity is approximately 1.0 fb^^. Assum- 
ing that the asymmetry A is due to asymmetric 5° 5° 
mixing and decay, we extract the CP-violation parameter 
of B'^ mixing and decay 0,|1I3: 
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■^12 (ri2) is the real (imaginary) part of the transi- 
tion matrix element of the Hamiltonian corresponding 
to {B^,B^) mixing and decay. Throughout this article 
we use the Particle Data Group 1] notation: B'^ ~ db, 
-Bg = sb. Ago is the dimuon charge asymmetry from 
direct-direct decays of B°B° (we define "direct decay" 
as 6 ^ H~X, and "sequential decay" as 5 ^ c ^ f^'^X). 
The dimuon charge asymmetry A in Eq. ^ excludes 
events with a muon from decay. Equation Q de- 
fines the factor /, to be obtained below, which accounts 
for other processes contributing to dimuon events. As a 
sensitive cross check, we also measure the mean mixing 
probability xo of B B hadrons, averaged over the mix 
of hadrons with a b quark. Finally we obtain the forward- 
backward asymmetry that quantifies the tendency of /x+ 
to go in the proton direction and /i~ to go in the anti- 
proton direction. 

The general case, with CP violation in both B'^ and 
systems, is considered in the last section of this arti- 
cle. In this general case, the dimuon charge asymmetry 
A has contributions from both 5° and B'^. Therefore, 
this measurement at the Fermilab Tevatron pp collider is 
complementary to similar measurements at B factories 
that are sensitive only to A^o , not A^o . 

The CP-violation parameter, defined in Eq. ((2Jl, is sen- 
sitive to several extensions of the standard model because 
new particles maycontribute to the box diagrams of M12 
[^Q. Reference Q concludes that "It is possible that the 
dilepton asymmetry could be one of the first indications 
of physics beyond the standard model" . 

The DO detector has an excellent muon system in Run 
II 0, with large (?y, 0) coverage, good scintillator-based 



triggering and cosmic ray rejection, low punch-through 
rate, and precision tracking. The muon is the particle 
with cleanest identification. The like-sign dimuon chan- 
nel is particularly clean: few processes contribute to it 
and fewer still contribute to an asymmetry. The DO de- 
tector is well suited for this precision measurement. 

The outline of the paper is as follows. The DO detec- 
tor is described in Section II. In Section III we consider 
the event selection. Physics and detector asymmetries 
are studied in Section IV. The processes contributing to 
the asymmetry A are presented in Section V, and their 
weights are summarized in Section VI. The breakdown of 
systematic uncertainties of A is discussed in Section VII. 
Cross-checks are listed in Section VIII. Final results are 
summarized in Section IX. 



II. THE DO DETECTOR 

The DO detector consists of a magnetic central-tracking 
system, comprised of a silicon microstrip tracker (SMT) 
and a central fiber tracker (CFT), both located within 
a 2 T superconducting solenoidal magnet 6]. The SMT 
has « 800, 000 individual strips, with typical pitch of 
50 — 80 /im, and a design optimized for tracking and ver- 
texing capability at pseudorapidities of \r]\ < 2.5. The 
system has a six-barrel longitudinal structure, each with 
a set of four layers arranged axially around the beam 
pipe, and interspersed with 16 radial disks. The CFT has 
eight thin coaxial barrels, each supporting two doublets 
of overlapping scintillating fibers of 0.835 mm diameter, 
one doublet being parallel to the collision axis, and the 
other alternating by ±3° relative to the axis. Light sig- 
nals are transferred via clear fibers to solid-state photon 
counters (VLPC) that have ~ 80% quantum efficiency. 

Central and forward preshower detectors located just 
outside of the superconducting coil (in front of the 
calorimetry) are constructed of several layers of extruded 
triangular scintillator strips that are read out using 
wavelength-shifting fibers and VLPCs. The next layer of 
detection involves three liquid- argon/uranium calorime- 
ters: a central section (CC) covering |?7| up to « 1.1, 
and two endcap calorimeters (EC) that extend coverage 
to 1 77! sa 4.2, all housed in separate cryostats In ad- 
dition to the preshower detectors, scintillators between 
the CC and EC cryostats provide sampling of developing 
showers at 1.1 < |?7j< 1.4. 

A muon system |5| is located beyond the calorimetry, 
and consists of a layer of tracking detectors and scintilla- 
tion trigger counters before 1.8 T iron toroids, followed by 
two similar layers after the toroids. Tracking at \ri\ < 1 
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relies on 10 cm wide drift tubes [7|, while 1 cm mini-drift 
tubes are used at 1 < [jyl < 2. 

A muon originating in a pp collision traverses the sili- 
con microstrip tracker and the scintillating fiber tracker 
in the 2 T solcnoidal magnetic field, the calorimeter, layer 
A of the muon spectrometer, the 1.8 T magnetized iron 
toroid, and layers B and C of the spectrometer. 

Luminosity is measured using plastic scintillator arrays 
located in front of the EC cryostats, covering 2.7 < < 
4.4. 

Trigger and data acquisition systems are designed to 
accommodate the high luminosities of Run II. Based on 
preliminary information from tracking, calorimetry, and 
muon systems, the output of the first level of the trigger 
is used to limit the rate for accepted events to ~ 2 kHz. 
At the next trigger stage, with more refined information, 
the rate is reduced further to « 1 kHz. These first two 
levels of triggering rely mainly on hardware and firmware. 
The third and final level of the trigger, with access to 
all the event information, uses software algorithms and a 
computing farm, and reduces the output rate to « 50 Hz, 
which is written to tape. 

The polarities of the toroid and solenoid magnetic 
fields are reversed roughly every two weeks so that the 
four solenoid-toroid polarity combinations are exposed 
to approximately the same integrated luminosity. This 
allows cancellation of first-order effects of the detector 
geometry. 



III. EVENT SELECTION 

Our standard cuts require global muons, i.e., local 
muon candidates (reconstructed from hits in layers A, 
B and C) with a matching central track (reconstructed 
from hits in the silicon and fiber trackers). To reduce 
punch-through of hadrons we only consider muons that 
traverse the iron toroid. To select muons that emerge 
from the toroid with momentum p > 0.2 GeV/c, we re- 
quire pt > 4.2 GeV/c or \pz\ > 6.4 GeV/c, where px is 
the momentum transverse to the beam measured by the 
central tracking system, and Pz is the component of the 
momentum in the direction of the proton beam. We re- 
quire at least two wire chamber hits in the A layer and at 
least three wire chamber hits in layers B or C. We require 
local and global track fits with good x^- To reduce cos- 
mic ray background, we require at least one scintillator 
hit associated with the muon to be within a time win- 
dow of ±5 ns with respect to the expected time. To re- 
duce muons from and tt^ decay, we require px > 3.0 
GeV/c. The track is required to have a distance of clos- 
est approach to the beam less than 0.3 cm. We use the 
full pseudorapidity range \t]\ < 2.2. We apply a cut of 
Pt < 15.0 GeV/c to reduce the number of muons recon- 
structed with wrong sign, and to reduce the background 
from and Z boson decay. This list completes the 
single muon cuts. 

The dimuon cuts are as follows. We require that both 



muon candidates pass within 2.0 cm of each other in the 
direction along the beam line at the point of closest ap- 
proach to the beam. To further reduce cosmic rays and 
repeated reconstructions of the same track (with differ- 
ent hits), we require the 3-dimensional opening angle be- 
tween the muons to be between 10° and 170°. We also 
require that the two muons have different A layer posi- 
tions (by at least 5 cm), different local momentum vectors 
(by at least 0.2 GeV/c), and different central track mo- 
mentum vectors (by at least 0.2 GeV/c). This completes 
the set of standard cuts. 

To avoid a bias due to mismatched central tracks 
(which are measured to be charge-asymmetric) we use 
the local muon charge, instead of the matching central 
track charge, to obtain the asymmetries. The positive 
charge asymmetry of central tracks is due to secondary 
particles that emerge from interactions with detector ma- 
terial. The measured charge asymmetry of central tracks 
with PT > 3.0 GeV/c is {N+ - N-)/{N+ + N') = 
0.0049 ±0.0005. 

The measurement of the dimuon charge asymmetry A 
is based on ratios of muon counts. To minimize the statis- 
tical error we use all recorded events regardless of trigger. 
If an event passes cuts, it is valid to accept it regardless 
of trigger, since an event with opposite muon charges can 
also be accepted by that trigger. 

Histograms oipT, rj, and (j) for standard cuts are shown 
in Figs, n and |21 Each plot superposes two histograms 
with opposite charge and opposite polarities (that can 
generally not be distinguished). 



IV. ASYMMETRIES 

In this Section we study detector and physics effects 
that may alter the charge asymmetry and therefore con- 
tribute to corrections and systematic uncertainties. The 
muon detector is shown schematically in Fig. |21 How can 
the detector introduce a charge asymmetry? It would 
have to have different acceptance x efficiency for posi- 
tive and negative muons, i.e., for tracks bending north 
and south in the magnetized iron toroid, see Fig. |31 
Such a difference may be due to an offset of the mean 
beam spot, to mechanical asymmetries, and to differ- 
ences in wire chamber and scintillator efficiencies. We 
find that the detector, operating with a given toroid and 
solenoid polarity, introduces an apparent dimuon charge 
asymmetry of approximately 0.006 in absolute value (to 
be discussed later in this Section). This detector effect 
changes sign when the toroid and solenoid polarities are 
reversed (since the exact same track that is called "posi- 
tive" with one polarity is called "negative" with the op- 
posite polarity). This effect can be seen by comparing the 
7] histograms in Figs. ^ and 12 Therefore, to cancel de- 
tector geometry effects to first order, we always consider 
data sets that have equal event counts for each toroid- 
solenoid magnet polarity (or weight the events appropri- 
ately) . We combine events with one solenoid polarity and 



6 




X10-* 



xlO 



2 4 6 8 10 12 14 
3 Pt [GeV/c] 




2 4 6 8 10 12 14 
Pt [GeV/c] 




50 
40 
30 
20 
10 




100 200 300 400 
3 [deg] 



50 
40 
30 
20 
10 




100 200 300 400 
<p [deg] 



FIG. 1: Distributions of single muon pt (a) and (d), r] (b) and (e), and (j) (c) and (f), for events with opposite toroid and solenoid 
polarities passing standard single muon and dimuon cuts. The charge times toroid polarity, af}, is indicated in the labels. Each 
plot superposes two histograms with opposite charge and opposite polarities. Positive (negative) charge corresponds to solid 
(dotted) lines (that can generally not be distinguished). Histograms with negative toroid polarity are scaled by the ratio of 
muon counts e = 1.0300, see Section IV. 



toroid polarity, with events with the opposite solenoid 
polarity and toroid polarity. The analysis is done sepa- 
rately for solenoid polarity equal to the toroid polarity, 
and solenoid polarity opposite to the toroid polarity. 

Let be the number of muons passing cuts with 
charge a = ±1, toroid polarity /3 = ±1, and 7 = +1 if 
77 > and 7 = — 1 if < 0. 

We model the physics and the detector as follows: 

ni'< = iiVe''(l+ayl)(l + Q7^/fc)(l + 7^d.t) 

X (1 + aP^Aro)(l + /37^/37)(l + olPA^p), (3) 

where e+ + = 1. The eight equations ^ define eight 
parameters in terms of the eight numbers . The pa- 
rameters are N , e^, and six asymmetries. Ne^^ is approx- 
imately equal to the number of muons passing cuts with 
toroid polarity (3. A is the dimuon charge asymmetry, 
A is the forward-backward asymmetry (that quantifies 
the tendency of /Lt+ to go in the proton direction and 
fjT to go in the anti-proton direction), Ad^t measures 



the north-south asymmetry of the detector ( "north" has 
77 < 0) , and Aro is the range-out asymmetry (that quan- 
tifies the change in acceptance and range-out of muon 
tracks that bend toward, or away from, the beam line, 
see Fig. O. A^p is a detector asymmetry between tracks 
bending north and tracks bending south. A and Aff, are 
physics asymmetries that we want to measure, and Adet, 
Aro and Aaf} are detector asymmetries. Apj is a second- 
order asymmetry that is different from zero only if A and 
Aro are different from zero. If the selection of events in- 
cludes single muon and dimuon cuts, we use capital A 
for the asymmetries. If only single muon cuts are re- 
quired, we use lower case a. In Tables HI and UTI we show 
the numbers n^'^ for our standard cuts. The measured 
asymmetries are presented in Table Hill 

We can understand the detector asymmetry Aap in 
more detail. Let A^ be the dimuon charge asymmetry of 
events with toroid polarity p. From Eq. ||3J) we obtain, to 
first order in the asymmetries, A^ = A+f3Aap. Therefore 
(3 Aaf} is the dimuon charge asymmetry due to detector 
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FIG. 2: Same as Fig. but for equal toroid and solenoid polarities. Here e = 0.9292. 
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FIG. 3: Schematic drawing of the magnetized iron toroids of 
the DO detector, and muon tracks related by toroid polarity 
reversal, CP conjugation and forward-backward reflection. 



geometry effects. It changes sign when the magnet po- 
larities are reversed. This is the detector geometry effect 
that is cancelled by taking the weighted average of A'^ 
and A~ . The magnitude of this effect is Aaf3 « —0.006, 



see Table Hm 

Solving equations ^ up to second-order terms in the 
asymmetries, A is obtained by taking the weighted aver- 
age + eA-): 



+ nl 



+ n_ ) 
^ A + AfbAdet,{4:) 

where e = /e~ . To the required accuracy, e is the ratio 
of the number of events passing cuts with toroid polarity 
(3 = +1 over the corresponding number with (3 = —1. 
The ratio e is determined by counting single muons for 
the single muon asymmetries, or by counting dimuons for 
the dimuon analysis. This procedure introduces no bias 
since we count all muons or dimuons regardless of the 
charges of each muon. The left-hand-side is the measured 
asymmetry, A is the corrected asymmetry, and —AfbAdet 
is the correction due to the forward-backward and detec- 
tor asymmetries. We do not apply this correction because 
it turns out to be negligible and compatible with zero, see 
Table ITlTl We use the correction to estimate the corre- 
sponding systematic uncertainty. We can understand the 
last term in Eq. Q): if positive (negative) muons prefer to 
go in the proton (antiproton) direction and the detector 
is north-south asymmetric, then we obtain an apparent 
charge asymmetry. 
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TABLE I: Numbers ng^ of muons passing standard single 
and dimuon cuts with charge a, toroid magnet polarity /3, 
and 77 < (7 = —1) or > (7 = +!)• There are two entries 
per event. The solenoid and toroid polarities are opposite. 
In total there are 154667 positive-positive, 154482 negative- 
negative, and 1075192 positive-negative dimuon events. No 
cuts on the dimuon mass are imposed at this stage. 



charge 


toroid 


-2.2 < 77 < 0.0 


0.0 < < 2.2 


a 


polarity (3 


7 = -l 


7 = +1 


+ 1 


+ 1 


367376 


335700 


-1 


+1 


348295 


353453 


+1 


-1 


337697 


343753 


-1 


-1 


356891 


325517 


TABLE 


IL Same as Table but for equal solenoid and 


toroid 


polarities. In 


total there are 


136422 positive- 


positive. 


, 136857 negative-negative, and 


944013 positive- 


negative dimuon events. 






charge 


toroid 


-2.2 < 77 < 0.0 


0.0 < < 2.2 


a 


polarity /3 


7 = -l 


7 = +1 


+1 


+1 


306594 


279508 


-1 


+1 


290270 


296264 


+1 


-1 


311153 


319602 


-1 


-1 


329285 


301908 



The forward-backward asymmetry A (up to second- 
order terms) is: 

(n|+ + nt' - nX~ - nl+) -I- e(n^+ + nZ~ - n" - 71I+) 
(n|+ + nt~ + + n++) + e(n^+ + nZ' + n:^_' + 7il + ) 

= Aft + AAdef{5) 



We have repeated the study of detector asymmetries 
for the central (I77I < 0.95) and forward (0.95 < \r]\ < 2.2) 
muon systems separately. 

Also shown in Table IIIII is the dimuon charge asym- 
metry for flavor creation Afc (defined as AiJ; > 90°) and 
flavor excitation Af^. (defined as Aip < 90°), where A-0 is 
the 3-dimensional angle between the two muons. Flavor 
creation corresponds to the b and b quarks in opposite 
jets, while flavor excitation corresponds generally to a 
bb pair produced in the hadronization of one parton, in- 
cluding gluon splitting. The accepted cross sections for 
flavor creation and flavor excitation are nearly equal (as 
indicated by their statistical errors in Table HTHl . In Ta- 
ble |lll| we also show the ratio R = {N++ + N — )/N+- 
of like-sign to opposite-sign dimuon events in either the 
mass window 5.0 to 8.7 GcV/c^ or 11.5 to 30 GeV/c^ for 
events passing a dimuon trigger. This ratio is used for 
mixing studies. These mass windows reduce backgrounds 
from same-side direct-sequential muon pairs (process P4 
in Table llV|) . and decays of J/^j's, T's and their reso- 
nances (part of process Pq), defined later. 



TABLE III: Asymmetries described in Section IV are shown 
for central (c) muons, I77I < 0.95, and forward (f) muons, 
0.95 < \ri\ < 2.2. All errors are statistical. The average for 
opposite and equal toroid and solenoid polarities (torpol-solpol 
= -1 and 1) is A = -0.0005 ± 0.0013 (all). 



torpol*solpol 


-1 


1 


e 


1.0300 


0.9292 


a (all) 


0.0001 ± 0.0005 


-0.0012 ±0.0005 


afb (aU) 


0.0021 ±0.0005 


0.0011 ±0.0005 


adet (aU) 


-0.0074 ± 0.0005 


-0.0045 ± 0.0005 


aro (aU) 


-0.0268 ± 0.0005 


-0.0298 ±0.0005 


a (c) 


-0.0014 ±0.0007 


-0.0035 ± 0.0008 


afb (c) 


0.0004 ± 0.0007 


-0.0010 ±0.0008 


adet (c) 


-0.0069 ± 0.0007 


-0.0068 ± 0.0008 


aro (c) 


-0.0867 ±0.0007 


-0.0891 ±0.0008 


a(f) 


0.0010 ± 0.0006 


0.0003 ± 0.0006 


aft (f) 


0.0033 ± 0.0006 


0.0025 ± 0.0006 


adet (f) 


-0.0078 ± 0.0006 


-0.0029 ±0.0006 


aro (f) 


0.0122 ±0.0006 


0.0089 ± 0.0006 


Afb (all) 


0.0006 ± 0.0006 


0.0001 ± 0.0007 


Adet (all) 


-0.0187 ±0.0006 


-0.0165 ±0.0007 


Aro (all) 


-0.0268 ±0.0006 


-0.0283 ±0.0007 


Acp (all) 


-0.0059 ± 0.0006 


-0.0069 ±0.0007 


Ai3^ (all) 


-0.0002 ± 0.0006 


-0.0015 ±0.0007 


R 


0.4603 ± 0.0013 


0.4610 ±0.0013 


Afc (all) 


-0.0008 ± 0.0026 


-0.0023 ±0.0027 


Afe (all) 


0.0019 ±0.0026 


-0.0008 ± 0.0028 


A (cc) 


0.0015 ± 0.0044 


-0.0021 ±0.0048 


A{S) 


0.0024 ± 0.0031 


0.0016 ±0.0033 


A (aU) 


0.0005 ±0.0018 


-0.0016 ±0.0019 



V. DIMUON PROCESSES 

In this Section we obtain the factor / defined in Eq. |5J). 
We consider the processes Pi - P13 listed in Table Hvl 
Pi is direct-direct bb decay. P2 is opposite-side direct- 
sequential decay. P3 is sequential-sequential decay. P4 
is same-side direct-sequential decay. P7 corresponds to 
cosmic ray muons that traverse the DO detector and are 
reconstructed twice: once upon entry and once upon 
exit. Pg corresponds to muons from decay, in co- 
incidence with a prompt muon from the collision. This 
process is discussed in detail in Sections VI and VII. 
Pg corresponds to a cosmic ray muon, in coincidence 
with a prompt muon. Pio corresponds to a hadron that 
traverses the calorimeter and iron toroid and is recon- 
structed as a muon, in coincidence with a prompt muon. 
Pii corresponds to a combinatoric background faking a 
muon, in coincidence with a prompt muon. Examples 
of "other" processes P12 are dimuons with the following 
parents: and tt^, P° and r*, P° and J/ifi, and 
T^, P" and J/4', B^ and r*, b and unrelated c. P13 are 
events that have one of the muons reconstructed with the 
wrong sign. Contributions from both muons coming from 
hadron misidentification or combinatoric background are 
negligible. 

Let Xd be the probability that a B^(bd) meson that 
decays to a flavor specific final state, mixes and decays 
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as a B°{db). Similarly, Xd is the probability that a 
meson mixes and decays as a B'^. Here we consider the 
possibility that Xd 7^ Xd (the general case with Xd 7^ Xd 
and Xs 7^ Xs is considered in Section IX). The probability 
that a b quark in the sample decays as a 6 is 



X — Id 1^0^ Xd + Is 1^0^ Xs 



(6) 



where fd and fs are the fractions of b quarks that 
hadronize to B^ or & and 5° or B^ respectively, and 
/3d, /3s and (/3) are the branching fractions for B", B^ 
and the 6-hadron admixture respectively decaying to 
with /X passing cuts. Similarly, 



r _P±- r Ps 

Id (^p^Xd + Js ^p^Xs 



(7) 



is the probability that a 6 in the sample decays as a 
b. From Ref. W\ we take fd = 0.397 ± 0.010, = 
0.107±0.011, Xdo = Uxd+Xd) = 0.186±0.004, and Xso = 
l(Xs + Xs) > 0.49883. We take (3d = Ps = (/3). To ab- 
breviate, we define xo = 5 (x + x), and ^ = 2xo (1 - XoJ- 
We assume CPT symmetry. The probability that a B (B) 
hadron that decays to a flavor specific final state, decays 
as a i? (B) is then 1 — xo- From Table ITVI we obtain the 
dimuon charge asymmetry A after correcting for asym- 
metric kaon decay, (i.e., after subtracting a term O.SaPg 
in the numerator), 

^ (X - X) [(1 - Xo)(A - ^3) + 0-25p'P^] ^ 



Ade 



and the factor / in Eq. ||2J): 



/ 



A 



den 



mXdo [(1 - Xo){Pi - P3) + 0.25p'P^] 



(9) 



where Aden = C{Pi + P3) + (1 - 0^2 + O.28P7 + 0.5P^ + 
0.78Pi3 and P;^ = Pg -I- Pg -t- Pio -t- Pn -f P12 . The fraction 
of prompt muons from b decay is p' = 0.6±0.15 i8|. 



VI. WEIGHTS P, OF DIMUON PROCESSES 

The weights P2 through P13, normalized to direct- 
direct bb decay Pi = 1, are summarized in Table 
The weights P2, P4, P5, Pe and P12 were obtained from 
Monte Carlo simulations (using the Pythia generator 
0) with full detector simulation (based on the Geant 
program 10]) and event reconstruction and selection. A 
cross-check for weight P2 is the measurement of the av- 
erage mixing probability of B hadrons to be described 
below. Weights P4, P5 and Pq do not contribute like- 
sign dimuons, and so do not enter into the measurement 
of the CP-violation parameter. Weight P3 was obtained 
from P3 « P2^/(4Pi). Weight P7 was obtained by two 
methods: (i) from the data of a cosmic ray run, and (ii) 
extrapolating the out-of-time muon background (as mea- 
sured by the scintillators) into the acceptance window of 



±5 ns. Weight Pg was obtained using the data of the 
cosmic ray run. Weight Pio was obtained by counting 
the number of tracks that had enough momentum to tra- 
verse the calorimeter and iron toroid, and multiplying by 
the probability exp (—14) that they do not interact (the 
calorimeter has ~ 7 nuclear interaction lengths, and the 
iron toroid has ~ 7 nuclear interaction lengths). Weight 
Pii was estimated by relaxing the number of required 
wire chamber hits (one less hit in layer A and/or 1 less 
hit in layers B or C). Weight P13 was estimated using the 
measured resolution of the local muon spectrometer. In 
our data set, passing standard single and dimuon cuts, 
we expect « 1 dimuon event from Z boson decay, and 
less than one event from prompt muons plus boson 
decay. 

Let us consider the weight Pg in some detail. This 
weight corresponds to prompt muons from b or c or s de- 
cay plus decay. This is an important background 
because kaon interactions are charge asymmetric, and 
dominate the systematic uncertainty of the measurement 
of the CP-violation parameter. The inelastic interaction 
length of K'^ in the calorimeter is greater than the in- 
elastic interaction length of . This difference is due 
to the existence of hyperons Y (strangeness —1 baryons: 
A, S, Y*). Reactions R-N Yn have no K+N ana- 
log. Therefore has more time to decay than K~. 
The result is a charge asymmetry from decay. The 
single muon charge asymmetry from decay is ob- 
tained from the inelastic cross sections for K^d and K'^d 
1] and the geometry and materials of the DO detector: 
a={n+- n-)/{n+ + n") = 0.026 ± 0.005. 

For Pg, we make two complementary estimates based 
on data. In the first, we measure the exclusive decay 
po ^ D*{2010)-p+i'f„ D*(2010)- D^t:^, D° 
K'^TT^ , and its charge conjugate. We apply standard 
single and dimuon cuts, and count events with and with- 
out a muon matching the kaon track. We subtract 
the background by two methods: using a side-band of 
m]j» - — niDO , or using the wrong relative sign of the muon 
from P° decay and the pion from D*~ decay. The result 
is Pg = 0.078 ± 0.013 (stat) ± 0.019 (syst). From studies 
with this exclusive decay, we learn that the global track 
X^-cut is not very effective in reducing decay kinks 
(for the high momentum muons passing cuts) . Therefore 
we must correct A for decay as discussed in Section 
VII. 

We use the following alternative procedure to esti- 
mate the background weight Pg from data. Instead of 
A'+ P^v^ we study 7r"'"7r~. For this estimate, 

we assume that the production and decay kinematics of 
Ar+ and A"^ are approximately the same. The branching 
fractions are similar. To account for the smaller decay 
length of AT^ compared to K^, we scale the volume in 
which the can decay (719 mm in radius and 1360 mm 
in half-length) by the fraction of lifetimes. 719 mm is the 
sum of the inner radius of the calorimeter, plus the trans- 
verse interaction length of K'^ in the calorimeter. We 
analyze single muon (instead of dimuon) data to account 
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TABLE IV: Processes contributing to dimuon events. Each row includes processes related by CP conjugation and B ^ B 
mixing. The weights are normalized to direct-direct 66 decay Pi = 1. ^ = 2xo (1 — Xo)- X = fdXd + faXe is the probability 
that 6 quarks decay as 6. x ~ fdXd + /sXs is the probability that 6 anti-quarks decay as 6. The fraction of prompt muons from 
6 decay is p' = 0.6 ± 0.15 8]. p = \p'(x ^ x)- = 0.026 ± 0.005 is the charge asymmetry of K"^ decay, see the text. CPT 
symmetry is assumed. For example, the number of direct-direct decays 66 p^ p^ X is oc Pix(l — Xo)- 



process 


weight 


N++ 




N+- 


6 6 ^ /i+ 


Pi = 1 


X(l - xo) 


(1 - Xo)X 


1-? 


6 p^ , b ^ c ^ p~ 


P2 


ia-c) 






b ^ c ^ p'^ , b ^ c ^ p~ 


Ps 


(1 - Xo)x 


X(l - xo) 




6 p~c /i"*" 


P4 








1 


c — > p'^ , c p~ 


Ps 








1 


Drell-Yan, J/ip, T 


Ps 








1 


dimuon cosmic rays 


Pr 


^ 0.14 


« 0.14 


« 0.72 


p + decay 


Ps 


0.25- (l + a + p) 


0.25- (1-a-p) 


0.5 


p + cosmic 


Pg 


R^O.25- (l-f/o) 


Si 0.25- (1-p) 


« 0.5 


p + punch-through 


Pio 


0.25- (1-fp) 


0.25- (1-p) 


0.5 


p + combinatoric 


Pll 


0.25 -(l-fp) 


0.25- (1-p) 


0.5 


other 


Pl2 


0.25- (l-fp) 


0.25- (1-p) 


0.5 


dimuon w. wrong sign 


Pl3 


0.39 


0.39 


0.22 



TABLE V: Weights of dimuon processes for standard cuts 
(obtained as described in the text). Note that 64% of dimuons 
are from direct-direct 66 decay. 



Pi 


= 1 




0.116 ±0.055 


Ps 


0.003 ±0.003 


Pi 


0.093 ± 0.049 


Ps 


0.070 ± 0.042 


Pe 


0.023 ± 0.023 


Pt 


0.003 ±0.003 


Ps 


0.078 ±0.023 


Pg 


0.0001 ±0.0001 


Pio 


0.001 ±0.001 


Pll 


0.0002 ± 0.0002 


Pl2 


0.163 ±0.066 


Pl3 


0.0005 ± 0.0005 



for correlations. We compare two histograms: one is pr 
of pions from if^'s decaying in the scaled-down volume, 
and the other histogram is px of the second (in order 
of decreasing pt) muon passing cuts. By this indirect 
method we obtain Pg « 0.041 ±0.010(stat)±0.041(syst). 

Within errors, the two measurements of Pg agree and 
we use the result from the first method. 

The systematic uncertainty of Pg, ±0.019, was taken 
as half the difference of 0.078 and 0.041. This value is 
reasonable in view of the variation of the measured Pg 
with different cuts and data sub-sets. 

From Monte Carlo simulations we obtain Pg = 0.047 ± 
0.034 (stat), consistent with the above. 



VII. SYSTEMATIC UNCERTAINTIES OF A 

We add in quadrature the following systematic uncer- 
tainties. A summary is presented in Table IVIl 

Detector effects. Before averaging over magnetic field 
polarities, the detector introduces a dimuon charge asym- 
metry of « 0.006 in absolute value, as discussed in 
Section IV. After averaging over magnetic field polari- 
ties (with appropriate weights), the uncertainty of the 
dimuon charge asymmetry due to detector geometry ef- 
fects is \Af}, ■ Adet\ (see Eq. Q). As a measure of this 
uncertainty, we have used the largest deviation from zero, 
\Afi,-Adet \ = 0.0049-0.030 = 0.00015, obtained after any 
of the 54 sets of dimuon cuts, for any of the three detector 
regions (central, forward, or all). 

Inaccuracy of e = /e~ . We have obtained the ratio 
e by counting dimuon events with toroid polarity /? = 1 
and dividing by the corresponding number for j3 — —\. 
This procedure introduces no bias since we count all 
dimuons regardless of the charges of each muon. We take 
Ae = 0.03 from the largest difference between any cuts. 
Multiplying by the detector dimuon charge asymmetry 
before averaging over magnet polarities, « 0.006 in ab- 
solute value, we obtain IS.A « 0.00018. 

Prompt /i + decay. The single muon charge asym- 
metry of kaon decay is a = 0.026 ± 0.005 as explained in 
Section VI. We take Pg = 0.078 ± 0.023 from Table |V| 
The corresponding correction to A, explained in Section 
V, is 5A = -Q.baPa/Aden = -0.5 x (0.026 ± 0.005) x 
(0.078 ± 0.023)/^de« = -0.0023 ± 0.0008 {Aden ~ 0.436 
is the denominator of Eq. ©). This uncertainty on 5 A 
is by far the dominating systematic uncertainty of the 
entire measurement. 

Dimuon cosmic rays. These are cosmic rays detected 
twice, once as they enter and once as they exit the DO de- 
tector. We take P-j < 0.007 from Table|V| From cuts that 
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TABLE VI: Systematic uncertainties of the dimuon charge 
asymmetry A for standard cuts. 



Source of error 


AA 






e = e+/6- 


0.00018 


prompt /i + decay 


0.00083 


dimuon cosmic rays 


0.00010 


prompt /i + cosmic fi 


0.00001 


wrong charge sign 


0.00018 


punch-through 


0.00001 


Total 


0.00089 


TABLE VII: Systematic uncertainties of /. 


Source of error 


A,f 


P2 


0.084 


Pr 


0.002 


PL 


0.056 


XdO 


0.012 


fd 


0.014 


XsO 


0.0002 


fs 


0.019 


P' 


0.008 


Pis 


0.0007 


Total 


0.105 



select cosmic rays, we obtain the apparent dimuon charge 
asymmetry A = —0.0095 ±0.0117. Then the correspond- 
ing uncertainty of A is < 0.007 x 0.28 x 0.021/Ade« = 
0.0001 (see Table ITvl 0.021 = | - 0.0095| + 0.0117). 

Prompt fi + single cosmic ray. We take Pg < 0.0002, 
see Table Ivl From cuts that select cosmic rays, we obtain 
an apparent single muon charge asymmetry a = 0.026 ± 
0.002. Then the uncertainty in ^ is < 0.0002 x 0.5 x 
0.028/Ade„ = 6 X 10-6 (see Table Hvl 0.028 = 0.026 + 
0.002). 

Wrong local muon sign. Pi3 < 0.001, see Table |V| 
Even if the asymmetry of wrong tracks is 0.1 (overes- 
timate), the corresponding uncertainty of A is small: 
< 0.001 X 0.78 X QA/Aden = 0.00018 (see Table Hvt. 

Punch-through. We take Pio < 0.002, see Table El 
The measured charge asymmetry of tracks with pT > 
3.0 GeV/c is 0.0049 ± 0.0005 due to showers on matter 
instead of antimatter. Then the error in A is < 0.002 x 
0.5 X 0.0054/^de« 1 X 10-^ (see Ta,h1e ITvl 0.0054 = 
0.0049 + 0.0005). 



VIII. OTHER CROSS-CHECKS 

The sign of central muons was cross-checked using cos- 
mic rays (which are charge asymmetric U). The sign of 
forward muons relative to the sign of central muons was 
cross-checked with J/tp's. 

The direct and reverse magnetic fields in the iron toroid 
were measured to be equal to within 0.1%. 



We find the dimuon charge asymmetry A stable 
(within statistical errors) for all recorded events or events 
passing a set of dimuon triggers, and across the differ- 
ent cuts (54 sets were studied), data subsets, opposite or 
equal toroid and solenoid polarities, central or forward 
muons, or flavor creation or flavor excitation events. 



IX. RESULTS 

We obtain A = 0.0005 ± 0.0018 (stat) for opposite 
solenoid and toroid polarities, and A — —0.0016 ± 
0.0019 (stat) for equal solenoid and toroid polarities (see 
last line in Table . Combining these measurements 
we obtain 



A = -0.0005 ±0.0013 (stat). 



(10) 



We add a correction 5 A = -0.0023 ± 0.0008 to A due to 
asymmetric decay (this effect is explained in Sections 
VI and VII). The uncertainty of this correction dominates 
the systematic uncertainties of the CP-violation parame- 
ter. The final corrected value of the dimuon charge asym- 
metry is 

A = -0.0028 ± 0.0013 (stat) ± 0.0009 (syst). (11) 

The breakdown of systematic uncertainties of A is pre- 
sented in Table ED 

From the dimuon charge asymmetry A we obtain 



A 



(12) 



-0.0023 ±0.0011 (stat) ± 0.0008 (syst). 



where 



/ = 0.814 ±0.105 (syst). 



(13) 



The breakdown of systematic uncertainties of /, calcu- 
lated from information provided in preceding sections, 
is listed in Table IVIII In comparison, the Particle 
Data Group average of 2004 [l| is 5R(eBo)/(l + |eBo|^) 
0.0005 ±0.0031. 

All preceding equations correspond to the case Xs — 
Xs- From ^ we obtain, for the general case Xd 7^ Xd and 



At 



fsXsO 
fdXdO 



A 



BO 



(14) 



We measure the ratio R of like-sign to opposite-sign 
dimuons. For this measurement we require the invari- 
ant mass of the two muons to be either in the window 
5.0 to 8.7 GeV/c2 or 11.5 to 30 GeV/c^. These mass 
cuts are designed to reduce backgrounds from same-side 
direct-sequential decay and backgrounds from ip and T 
meson decays to allow a measurement oi B B mix- 
ing. For this reason we set P4 = Pg = for the mix- 
ing analysis. We also require dimuon triggers from a 
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TABLE VIII: Systematic uncertainties of ^ = 2xo(l — Xo)- 

Pl = P4 + P5 + Pe. 

Source of error 

R 0.0074 

P2 0.028 

Pi 0.015 

Pr 0.0001 

Pi 0.015 

Pi3 0.0002 

Total 0.036 



A fb is defined in Section IV. The systematic uncertainty 
is |A||Adet| < 0.0044 X 0.030. As indicated earlier, the 
fraction of muons from W decay in this sample is negli- 
gible. 

In conclusion, the results 1)111) through 116(1 . are consis- 
tent with the standard model 1], and constrain some of 
its extensions The general result H14|l complements 

measurements at B factories, which are sensitive only to 
Ago, not Abo. 



list that excludes triggers requiring opposite sign muons. 
We obtain R = 0.461 ± 0.001 (stat) ± 0.010 (syst), and 
^ = 0.229 ± 0.001 (stat) ± 0.036 (syst). The breakdown 
of systematic uncertainties, calculated from information 
provided in preceding sections, is shown in Table IVIIII 
The final result for the mixing probability, averaged over 
the mix of hadrons with a b quark, is 

Xo = 0.132 ± 0.001 (stat) ± 0.024 (syst). (15) 

The agreement with the world average, 0.127± 0.006 0, 
is a sensitive test of P2 and /, since the largest systematic 
uncertainty of ^ and / are due to the same weight P2 . 

Finally, we measure the tendency of /i+ (/i^) to go 
in the proton (antiproton) direction. We obtain the 
forward-backward asymmetry for events passing stan- 
dard single and dimuon cuts: 

Afb = 0.0004 ± 0.0005 (stat) ± 0.0002 (syst). (16) 
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